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ABSTRACT 

We examine the expected X-ray polarization properties of neutron-star X-ray 
sources of various types, e.g., accretion and rotation powered pulsars, magnetars, 
and low-mass X-ray binaries. We summarize the model calculations leading to 
these expected properties. We describe how a comparison of these with their ob- 
served properties, as inferred from GEMS data, will probe the essential dynam- 
ical, electromagnetic, plasma, and emission processes in neutron-star binaries, 
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discriminate between models of these processes, and constrain model parame- 
ters. An exciting goal is the first observational demonstration in this context 
of the existence of vacuum resonance, a fundamental quantum electrodynamical 
phenomenon first described in the 1930s. 



1. Introduction 

A major class of X-ray sources that the Gravity and Extreme Magnetism SMEX (GEMS) 
would be studying is neutron stars in various astrophysical situations, operating as accre- 
tion powered pulsars (APPs), rotation powered pulsars (RPPs), magnetars, accretion pow- 
ered millis econd X-ray pulsars (AMXPs), and brigh t, non-pulsing low- mass X-ray binaries 



(LMXBs) fIShapiro fc Teukolskvl Il983l : iGhoshI l2007f ). APPs, AMXPs and LMXBs are bi- 



nary systems where a neutron star accretes matter from its companion, and conversion of 
the gravitational binding energy of the accreted matter powers X-ray emission. RPPs are 
thought to generate X-rays by one or more of the mechanisms summarized in Secj3l the 
ultimate source of energy for this emission being that of the neutron star's rotation. The 
majority of RPPs are single, relatively young neutron stars, but there is a population of 
RPPs in binaries, consisting mostly of old, recycled neutron stars with low-mass compan- 
ions, and also a small number of young neutron stars with massive companions. Magnetars 
are believed to generate X-rays by conversion of the magnetic energy available from the su- 
perstrong magnetic fields of the neutron stars in them. These objects are usually thought to 
be single, although there have been occasional suggestions of magnetars in binaries. All of 
the above classes of neutron-star X-ray sources exhibit periodic pulses at the rotation period 
of the neutron star, except for LMXBs, for which these pulses are believed to be unobservable 
or absent for reasons summarized in SecJSl although the neutron stars in these sources are 
thought to be rotating fast as a result of the recycling process which they are undergoing. 
However, many LMXBs undergo X-ray bursts, and during the bursts from some of them, 
periodic oscillations have been observed at high frequencies, which are believed to be the 
spin frequencies (or simple multiples/sub multiples thereof) of the underlying neutron stars. 
A classic timing property exhibited by many LMXBs is quasiperiodic oscillations (QPO) at 
low, medium, and high frequencies, to which we return in Sec J6.2[ Correlating the above 
timing behaviors of the X-ray fluxes from these classes of neutron-star sources with those 
of their X-ray polarization properties (amount /degree of polarization and its position angle, 
or, alternatively, the Stokes parameters) will be a principal diagnostic available to GEMS 
for addressing the scientific questions listed in SecH 



The magnetic field of the neutron star is the most important parameter for categorizing 
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the X-ray polarization properties at the source [i.e., the environs of the neutron star), as 
expected, since the primary cause for the emission polarized X-rays from APPs, RPPs and 
magnetars is this magnetic field — directly or indirectly — as explained below. Among 
the above classes of neutron-star sources, APPs and RPPs are universally believed to have 
the "standard" surface field-strengths of ~ 10^^ G canonical for neutron stars, with a rough 
spread of about one order of magnitude around this value. Magnetars are thought to have 
superstrong fields ~ 10^^ — 10^^ G by this standard, whence the name. By contrast, LMXBs 
and AMXPs are thought to have fields ~ 10^ — 10^ G, low by the above standards. Such 
low fields are believed to have been produced by a reduction of the standard neutron-star 
field stre ngths in the process o f recycling which occurs during the long accretion phases in 



LMXBs fivan den Heuvellll992[ ) 



For APPs and magnetars, it is their magnetic fields that dominate both (a) the produc- 
tion processes of polarized X-rays at or near the neutron-star surface, and, (b) the subsequent 
propagation of these polarized X-rays through the magnetosphere surrounding the neutron 
star (see Secsj2] and Hj). This happens because the magnetic field strongly influences both 
(a) the differential opacity between the ordin ary (O) and extraordina. ry (E) polarization 



modes, and, (b) the scattering cross-sections (IMeszaros fc Nagell ll985al Jbl). For RPPs, the 
X-ray emission process in the magnetosphere is dominated by the local magnetic field, and so 
therefore is the polarization of the resultant radiation (Secj3]). For LMXBs and AMXPs, by 
contrast, the magnetic field is so low that its effects are thought to be negligible in the con- 
text of X-ray polarization, so that the latter is described in terms of non-magnetic Compton 
scattering (see SecsE]and[H]). 

In addition to its magnetic field, the other essential property which characterizes a 
neutron star is its spin period, so that it is most useful to display neutron stars in a period 
- magnetic field (P — B) diagram, as shown in FigJU The different classes of neutron-star 
systems introduced above occupy characteristically different areas in this diagram, as shown. 
The dominant X-ray polarization-producing mechanisms are also characteristically different. 
For APPs and magnetars, it is the strong differential opacity between O and E modes, while 
for RPPs, it is the inherent strongly-polarized nature of the synchrotron/curvature radiation 
emitted from their magnetospheres, and in both cases the polarization degree can be very 
high, upto ~ 80%. For LMXBs and AMXPs, it is the characteristically low polarization 
(~ 10% — 15%) signature given to the radiation by non-magnetic Compton scattering in the 
accreting plasma around the neutron star. 

In the following sections, we summarize the current state of knowledge of the expected 
X-ray polarization properties of the above classes of neutron-star sources, and GEMS strate- 
gies for observing and utilizing them for diagnostics of emission geometry, magnetic field 
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structure, and fundamental quantum electrodynamic (QED) effects. In Secj2l we consider 
APPs. RPPs are considered in Secj3l and magnetars are dealt with in SecJH Low magnetic 
field systems are considered in later sections, SecE] dealing with AMXPs and SecE] describ- 
ing LMXBs, including both burst oscillations and QPO sources. Finally, we summarize our 
conclusions in SecJU discussing the key scientific questions to be addressed by GEMS. 



2. Accretion Powered Pulsars (APPs) 



Plasma accreting onto the neutron star in an APP is channeled by the stell ar magnetic 



field t o the magnetic poles and s o forms accretion columns above these poles ( iLamb et al. 



19731 : iDavidson fc Ostrikerlll973l ). X-rays generated at the bases of the accretion columns 
are transported through the highly magnetized plasma there, and the transport properties 
{e.g., the scattering and absorption cross-sections) depend strongly on the polarization of 
the X-rays, being very different for the ordinary (O) mode and the extraordinary (E) mode. 
When vacuum polarization effects (see below) are negligible, the O mode has its electric 
vector (which is the conventional definition of the direction of polarization) lying in the 
plane defined by the external (z.e., stellar) magnetic field B and the wave vector k, while 
the E mode has its direction of polarization perpendicular to the (k , B) plane. When 
vacuum polarizatio n effects are significant, the mode characteristics are more complicated 
(jPavlov et a/.lll980l ). This strong disparity between the transport of O and E modes leads to 
a strong polarization of the emergent radiation. A simple, rule-of-thumb argument widely 
used to illustrate the point is that, at photon energies E = hv far below the cyclotron 
energy E^ = h{eB/mec) ~ II.6-B12 keV (here, B12 is B in units of 10^^ G), the opacity ke 
of the E mode is drastically reduce d compared to that of the opacity kq of fhe O mode, 
scahng roughly as k^; ~ {E / EcYkq (ILai fc Hdl2003l ). Consequently, the emergent radiation 
is dominated by the E mode, which comes from the deeper and hotter layers of plasma, 
and so is strongly polarized in the direction of the E mode. We shall see below that the 
actual situation is more complicated, but that the general thrust of the above argument is 
qualitatively correct. 

Emission of X-ray pulses from APPs is conceptually visualized in terms of rotating X- 
ray beams which can be characterized as either pencil- or /an-shaped in two useful limits. 
In the former, the opacity along the accretion column is sufficiently low that the beam 
emerges preferentially along the column, in a pencil shape. In the latter, the opacity along 
the column is so high that the bea m emerges preferentially through the sides of the column, 
i n a fan shape (Lamb a/.l Il973l ). Pulse profiles for these beam g eometries are d i fferen t 
(IMeszaros fc Nagell Il985bl ). and so are the polarization properties (IMeszaros et a/.lll988l ). 
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We follow here the description of these properties pioneered by Meszaros and co-authors on 
the basis of their extensive numerical treatr aent of radiative transfer through the strongly 
mgnetized plasma in the accretion columns ( iMeszaros fc Nagellll985al ). 



In the Meszaros et al. approach, the pencil beam is modeled as emission from the face of a 
plasma slab lying on the stellar magnetic pole (the stellar magnetic field being perpendicular 
to the slab surface), while the fan beam is modeled as emission from the sides of a plasma 
cylinder sticking out of the stellar surface around the magnetic pole (the stellar magnetic field 
being along the cylinder axis). Transfer of radiation is handled in either geometry through 
the Feautrier formalism (jFeautrierlll964l ). using scattering and absorption cross-sections for 
O and E modes, and a redistribution matrix that keeps track of the scattering of photons 
from one angle to another, from one energy to another, and from one polarization mode to 
the other (IMeszaros fc Nagellll985al ). 



In a strong magnetic field, not only is the plasma birefringent, but so also is the vac- 
uum, the latt er being a fundamental qu antum electrodynamical (QED) effect first pointed 
out long ago (IHeisenberg fc Euleii Il936l ). These two effects generally work "against" each 
other in a sense, and so "compensate" each other at a point called vacuum resonance 
( Meszaros fc Ventura 1979 ). where the degree of linear polarization vanished. This vac- 
uum resonance energy depends on both the magnetic field strength B and the electron 
density ne, and is given by ^ 133^2 n'],''22 keV (here, ne.22 is in units of 10^^ cm~^, and 
B12 as before). For photon energies below E^, the normal modes of the plasma dominate 
in the radiation-transfer process, while those of the vacuum dominate at E > Ey. Cross- 
sections for scattering and absorption for both E and O modes sh ow strong features both a t 
the cyclotron energy Ec and at the vacuum resonance energy Ey (IMeszaros fc Nagellll985al ). 
so that the polarization properties of the emergent X-rays from APPs go through strong 
changes when the X-ray photon energy passes through either of these resonances, providing 
valuable diagnostics of magnetic fields and plasma conditions. For canonical field strengths 
of APPs (see above), it is clear that Ec would be outside the energy range of GEMS op- 
eration, so that our focus here will necessarily be on the expected changes in the essential 
polarization parameters at Ey. 

The GEMS energy-band of ~ 1 — 10 keV defines an allowed region in the Magnetic 
Field - Luminosity {B — L) plane of APPs, so that, within this region, Ey lies inside the 
GEMS band, and therefore valuable diagnsotics at Ey would be possible with GEMS. A 



^Actually, the directions of the dominant linear polarizations generated by plasma and vacuum birefrin- 
gence are perpendicular to each other. At vacuum resonance, their strengths are equal, so that the resultant 
polarization is circular, not linear. 
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rough sketch of this allowed region is shown in Figj2l fo r canonical properties of accreting 



neutron stars and accretion geometry (iGhosh et a/.ll2013l ) 



Model X-ray pulse profiles of APPs are shown in Figs|3]and|ll together with the (linear) 
X-ray polarization degree Pl and the polarization angle x as functions of the pulse phase 0, 
as obtained from the numerical scheme of Meszaros et ai, referred to henceforth as the M 
scheme. For the particular neutron-star and accreting-plasma parameters used in the results 
displayed in FigsE] and SJ the cyclotron resonance energy is Ec ~ 38 keV, and the vacuum 
resonance energy is fa 22 keV. (Note that we have shown here for illustration the original 
M scheme results obtained for the parameters used by those authors. Recent estimates of 
these parameters for well-known, bright APPs i ndicate that -E,, is expected to be in the 



GEMS band for a substantial number of them (iGhosh et a/.ll2013[ ).) We shall henceforth 
call libraries of such model profiles of total X-ray intensity (/), P^, and x as functions of 
essential system parameters {e.g., the inclination angle and the observer angle) as atlases. 
Comparison of such atlases with GEMS data on APPs will be a major diagnostic probe in 
our work, and we hereby name this probe pulse phase polarimetry, which will indeed be a 
most useful probe for X-ray pulsars of any kind, as we shall see in subsequent sections. 

A key signature of the beam geometry (pencil or fan) is the phase correlation between 
the pulse profile (which we can also call the intensity profile or /-profile) and the P^-profile, 
which can be seen clearly in Figsj3]and|H For pencil beams, the maximum in Pl is generally 
in phase with the pulse maximum for photon energies below E^, but generally out of phase 
for photon energies above Ey. For fan beams, the situation is exactly the opposite. In a 
similar vein, the phase correlation between the pulse profile and the x-profile is also a good 
diagnostic. For pencil beams, x generally goes from positive to negative values at the pulse 
maximum, and from negative to positive values at the pulse minimu Again, for fan 
beams, the situation is exactly the opposite. 

It follows from the above that a good pictorial representation of the I — Pl diagnostic 
would be a direct I — Pl plot, which can be obtained from the /- and P^-profiles, treating 
the phase as a parameter. This I — Pl plot would be a closed, ellipse-like, but more 
complicated curve, in analog y with the familiar Lissajous figure which would occur if the 



above profiles were sinusoids (IGhosh et a/.ll2013l ). The major axis of this figure would have 



a positive slope below E^, but a negative slope above this energy for pencil beams, so that 
the figure would turn dramatically by about 90 degrees as the observation energy passes 



^ We exclude here discontinuous jumps x by 180 degrees, which occur when the line of sight points directly 
down at the magnetic pole at some pulse phase, and there is a discontinuous jump in the angle between 
sky-projected stellar magnetic field vector and the sky-projected rotation axis. 
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through Ey. For fan b eams, the situation would be the opposite, and the turn would be in 
the opposite direction (iGhosh et a/.ll2013l ). Similar considerations apply to the I — x plot. 



The slab or cylinder shaped em ission regions envisaged in the above M scheme are typ- 
ically of radial extents ~ 1 — 10^ m (IMeszaros fc Nagellll985al Jbl). i.e., small compared to the 
neutron-star radius, and quite tiny compared to the typical extent of an APP magnetosphere 
(~ 10^ — 10^ m). Polarized X-rays produced by these regions propagate subsequently through 
the magnetospheres of APPs, which lead to possible modifications of their properties. There 
are two major propagation effects whose polarization signatures on X-rays from APPs need 
to be quantified fully. 

The first effect is the resonant Compton scattering of an outgoing photon generated by 
the M scheme from those points in the accretion column where the photon's frequency equals 
the local cyclotron frequency in the rest frame of the infalling plasma. In the rest frame of 
the neutron star, then, the corresponding condition in terms of the energies defined earlier 
can be written as 

E=-.^. (1) 



7(1 - M ■ 

Here, /i = k.B is the cosine of the angle between the photon propagation (k) and magnetic 
field (B) directions, /3 = f/c is the dimensionless electron velocity, and 7 = (1 — is the 
Lorentz factor. There is a limiting surface of last scattering, or escape surface, outside which 
it is not possible to satisfy the condition given by Eg. ([I]). The location o f this escape surface 



is given by the condition (Ec/E)"^ + /j,"^ = 1 (IFernandez &: David I2OIII ). Its radius r^s has 



been calculated for magnetars (see SecH]), and the values for APPs with characteristically 
lower magnetic fields (see Sec{T]) are correspondingly lower, being given for a purely dipolar 
neutron-star magnetic field roughly as 



R 



2.3 



NS 



E \ 



1 1/3 



IkeV 



) B,2 (1 - /i') 



2\-l/2 



(2) 



The second effect is that of the polarization freezing radius Vpi. As the X-ray photons 
propagate outward in the magnetosphere, normal modes do not mix significantly at first: 
this propagation is called adiabatic, as each mode follows its own evolution in the changing 
magnetic field without being affected by the other one. At a sufficiently large radius Vpi and 
small magnetic field strength, however, this adiabaticity breaks down and the modes mix, 
thereby fixing the polarization degree and direction for all further outward propagation, rpi 
has been calculated for magneta rs (see SecHl), where the e ntire magnetosphere is dominated 
by vacuum polarization effects (IFernandez fc DavisI 120111 ). This can be seen by calculating 
the radius of the vacuum resonance point (see above), and remembering that vacuum 
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polarization effects dominate for r < r.^, wliile plasma polarization effec ts dominate for r > 



For magentars, ~ SOOORjys, where R^s is the neutron-star radius (IFernandez fc Davis 



201ll ). For APPs, on the other hand, the magnetosphere is almost entirely d ominated by 



plasm a polarization effects, since the vacuum resonance radius is given by (iGhosh et al. 
2013h 

^ YsiTRri^] ' i^s-rv (3) 



Rns 

It is clear from Eq. ( 



0.91 



(. 



19/49 



VlKevy "° V^0. 

that, except at high photon energies, high magnetic fields, and 
very low luminosities, the entire magnetosphere is dominated by plasma polarization effects. 
Thus, a calculation of Tpi for APPs involves only the properties of plasna a polarization modes. 



and so is entirely different from that for magentars (iGhosh et al. 



2013|). 



A third, interesting, fundamental effect is that of general relativity on photon propaga- 
tion near the neutron-star surface, but it appears that its strength is not large enough for 
neutron-star radii given by modern equations of state (EOS) to warrant a detailed inclusion 
into a polarization computational scheme at this stage. Gravitational bending of photon 
trajectories ca n lead to interesting effect s on the yulse profiles of A PPs which were studied 
in the 1980s feieffert fc MeszarosI Egssi : iMeszaros fc RieffertI Egssh . Briefly, a fraction of 
the photons emitted near the stellar surface can be so bent as to propagate backward, so 
that this backward flux would partly (a) be bl ocked by the stellar surfac e, and partly (b) 
enhance the flux from the other magnetic pole (IRieffert &: MeszarosI Il988[ ). This would oc- 
cur for both pencil and fan beams in the M scheme introduced above. Consequent changes 
in the p ulse proflle, particularly in the degree of modulation, were studied by the above 
authors ( iMeszaros fc Rieffertlll988[ ). As expected, the strength of the GR effects increases 
with increasing compactness of the neutron star, i.e., decreasing values of the parameter 
a = Rns/Rs, where Rs = 2GMns/c^ is the Schwarzschild radius. For M^s = 1.4M0, 
Rg 4.2 km. The above authors studied the parameter space 1.6 < a < 4.0 (which corre- 
spond to Rns between 6.7 km and 16.8 km) and showed that GR effects became important 
for a < 2.0, which corresponds to R^s < 8.4 km. Those neutron-star EOS which are 
currently considered viable give values of Rns in a narrow range around ~ 12 km (corre- 
sponding to a ~ 3), where GR effects would b e minor, estimated t o be ~ 10% or less, for 
both pulse profiles and polarization properties (IMeszaros et a/.lll988l ). Accordingly, we shall 
not consider these effects further in this paper. 
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3. Rotation Powered Pulsars (RPPs) 

Radio polarization measurements have been a valuable probe of the emission geometry 
in RPPs, and the only measurement {i.e., not an upper limit) of X-ray polarization from 
neutron stars so far h as been from the well-known RPP, the Crab pulsar, and its nebula 



(IWeisskopf et a/.lll976l ). However, the emission sites and mechanisms for X-ray and radio 
emission are believed to be entirely different because the X-ray and radio pulse profiles are 
quite different in general. There are three main types of models for X-ray (and other high en- 
ergy) emission from RPPs, namely, (a) polar cap models (see Daugherty & Harding 1996 and 
references therein), where the em ission occurs withi n a few stellar radii of the neutron-star 



surface, (b) the outer gap models (ICheng et a/.lll986l ). where the emissio n occurs in the oute r 



magnetosphere near the light cylinder, and, (c) striped wind models (iPetri fc Kirkll2005l ) 
where the emission occurs in the pulsar wind outside the light cylinder. In the slot gap model 
a recently explored variation of polar cap-type models, the emission o ccurs at high altitudes 



along the last open magnetic field line, in the outer magnetosphere ( iMuslimov fc Harding 



20031 ). 



Polar-cap models assume that particles begin accelerating near the neutron-star surface, 
and emit high-energy radiation at radii comparable to -Rats' through curvature radiation 
and/or inverse Compton- induced pair cascade in the strong magnetic fields there. The slot- 
gap model is somewhat similar, except that the curvature, synchrotron, and inverse Compton 
components of the radiation originate in the outer magnetosphere. Outer-gap models assume 
that the acceleration occurs in vacuum gaps that form between the null-charge surfaces and 
the light cylinder in the outer magnetosphere, and that the observed high-energy radiation 
is curvature radiation and photon-photon pair-production induced cascades from these sites. 
Striped-wind models envisage the pulsed emission as originating in the pulsar wind beyond 
the light cylinder, due to toroidal magnetic fields of alternating sign in this wind, which 
can form current sheets for pulsars with high inclination, and the reconnection of these 
magnetic fields can convert magnetic energy into particle energy, leading to sharply (double) 
peaked pulse profiles if the dissipation region is small enough. Both slot-gap and outer-gap 
models have caustics in their emission pattern, i.e., photons emitted over a broad range of 
altitudes tend to pile up at roughly the same phase, due to an almost complete cancellation 
between (a) the phase delay in photon emission with increasing altitude along a trailing field 
line, and, (b) special-relativity effects like the aberration of photon-emission directions and 
the time-of-fiight effects caused by the finite speed of light. In the two-pole- caustic (TPC) 
configuration, strong caustics form on field lines trailing both polar caps, leading to light 
curves with widely-separated sharp double peaks. 

In the light of extensive studies of high-energy emission from RPPs with Fermi over the 
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last several years, it is now clear that observations strongly disfavor polar-cap type models 
where this emission occurs near the neutron-star surface, and strongly support those types 
of models in which this emission takes place in the outer magnetosphere or beyond the 
light cylinder (Harding 2010 and references therein). Accordingly, we shall consider only 
the latter types here for X-ray polarization studies, giving an occasional reference to results 
for polar-cap type models to show the difference. Expected pulse-phase variation profiles 
of the X-ray polarization have to be calculated numerically for these models, since simple, 
empirical models like the rotating vector model, so useful for lower-altitude radio emission 
(see, e.g., Ghosh 2007 and references therein), are no longer relevant here. Results from 
such calculations for the slot-gap and outer-gap models are shown in FigJSl with that for 
the polar-cap model also shown for reference. The latter shows an S-shaped swing in the 
polarization angle, reminiscent of the rotating vector model, but irrelevant here for the above 
reasons and also in strong disagreement with optical polarization measurements of the Crab 
pulsar. The polarization signatures of the TPC/slot-gap and outer-gap models are very 
different indeed. In particular, for the TPC/slot-gap models, the polarization angle shows 
very rapid sweeps through the main peak and the interpulse, and a strong depolarization 
associated with each peak. These features a. re actually observed in the optical polarization 



profile of the Crab pulsar (ISmith et a/.lll988[ l. 



As the different outer-magnetosphere models produce high energy RPP pulse profiles 
which are rather similar, it would be difficult to discriminate between them on this count. 
However, since different models of this type yield quite different polarization profiles at these 
energies because of their different magnetic field structures in the outer magnetosphere, pulse- 
phase polarimetry in the X-rays would serve as an excellent discriminator between models, 
and this would be a major direction of RPP study with GEMS. To this end, we shall have 
atlases of the profiles of J, Pl, and x corresponding to the above models, each such atlas 
detailing the profiles for a particular model for a grid of values for the essential parameters, 
e.g., the inclination angle and the observer angle. The idea is exactly analogous to what we 
described in Secj2] for APPs, and examples of such RPP atlases are shown in FigslG] and [3 



4. Magnetars 

Magnetars are believed to generate X-rays by conversion of the magnetic energy stored 
in the twisted magnetospheres of neutron stars with superstrong magnetic fields (see SecH]). 
This energy release near the stellar surface is thought to produce "hot spots", which emit 
thermal X-rays through the very thin (typically 1 — 10 cm thick) neutron-star atmosphere. 
The polarization properties of this primary radiation are determined by the interplay be- 
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tween the radii of the photospheres of O and E modes on the one hand, and the vacuum 
resonance radius r„ on the other (see Lai & Ho 2003 and references therein). At magnetar 
field strengths, the vacuum resonance hes between the O- and E-mode photospheres (the 
latter always lying deeper because of the lower opacity of the E-mode), and the emergent 
radiation from the atmosphere is dominated by the E mode at all photon energies of in- 
terest for GEMS observations. Expected pulse-phase var iations in the polarization of this 
radiation, as given by model calculations ( iLai fc Hoi 120031 ) . are shown in FigJHl 



Propagation of the above polarized X-rays through magnetar magnetospheres has two 
major effects on the polarization, which have been already introduced in SeclH viz., resonant 
Compton scattering and polarization freezing. The first is described by the limiting surface 
of last scattering or the escape surface, outside which it is not possible to satisfy the condition 
given by Eq. ([T]). The location of this esc ape surface is given by the condition {Ed E)"^ + ^'^ = 
1 and its radius res has been calc ulated (iFernandez fc David |20 111 ) for the standard twisted- 
dipole magnetospheric structure ( iThompson et a/.ll2002l ) of magnetars to be 



R 



12 



NS 



E 



IkeVy 



-1/2 



l/2.i 



(4) 



Here, B14 is the polar magnetic field of the magnetar in units of 10^^ G, ^ is a geometrical 
factor describing the angular dependence of the magnetic field, and t he specific result given 
by Eq.(|l]) is for a twist parameter of unity ( iFernandez fc David 1201 ll ). 



The second effect, viz., polarization freezing, is described in terms of a polarization 
freezing radius Vpi, the physics of which has been introduced in SecJJl Magnetar magneto- 
spheres are completely dominated by vacuum polarization effects, as seen by calcula ting the 
vacuum resonance point radius (see Secj2]) r„ ^ 3000-RAr5 (iFernandez fc David l201l[ l. Using 
the properties of the vacuum modes, Vpi for magnetars has been calculated as 



R 



146 



NS 



-1 1/4.76 



(5) 



Here C is a dimensionless function of order unity involved in calculating the scale length on 
which B changes in the magnetosphere. 

The above separation by roughly one order of magnitude between -Rats', r^s, and Vpi 
for magnetars makes the calculation of propagation effects relatively straightforward for 
them, since resonant scatt ering is basically decoupled from subsequent polarization freezing 
(IFernandez fc David I2OIII ). and both are, in a sense, decoupled from the processes in the 
neutron star's atmosphere. Such calculations have been done for the propagation of an 
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assumed completely polarized radiation from the neutron-star atmosphere through a model 
magnetar magnetosphere, and the results for the final emergent radiation are summarized in 
Figini The propagation effects tend to reduce the amount of polarization, but polarization 
degrees are still relatively high (~ 40 — 80%) at lower (~ 2 — 4 keV) energies, and particularly 
so around the pulse minima. The calculated pulse-phase variation profiles of the X-ray 
polarization of magnetars are generally similar to those of pencil-beam APPs (described 
by the slab model in the M scheme: see Secj2]), as expected, since hot-spot emission from 
neutron-star surfaces in magnetars would have a pencil-beam character. 



5. Accretion Powered Millisecond X-ray Pulsars (AMXPs) 



In this and the next section, we consider systems containing accreting neutron stars 
with low magnetic fields ~ 10^ — 10^ G, where the magnetic field is assumed to have a 
negligible effect on radiation transport, and X-ray polarization properties of such systems 
are described in terms of Compton scattering in a plasma configuration with a geometry 
appropriate for the system. We consider AMXPs first, of which there are about thirteen 
currently known. These systems are believed to have been produced by the recycling of 
neutron stars in low-mass X-ray binaries (LMXBs: see next section), where accretion had 
both (a) spun up the neutron stars to millisecond periods and (b) reduced their magnetic 
fields to the above levels. 

The emission sites of AMXPs are modeled as "hot spots" on the neutron star surface, in 
analogy with APPs. Thermal radiation from these hot spots undergoes Compton scattering 
in the plasma behind the accretion shock, which is modeled as a plane-parallel slab of optical 
depth of ord er unity (r ~ 1). The basic polarization signatures of Compton scattering are 
well-known (I Chandr asekharl Il960l : iLoskutov Sz Sobolevi Il982l : iNagirner fc Poutanenlll993h . 



For an optically thick (r >> 1) slab the direction of the linear polarization is pref- 
erentially in the plane of the slab, and for the classic problem of a semi-infinite slab in 
the Thompson scatterin g limit, the maximuna degree of polarization possible is the Chan- 
drasekhar limit ~ 12% (j Chandr asekharl Il960l ) . Detailed calculation s for r = 1 have been 
perfo rmed, which yield typical polarization degrees ~ 10% — 15% (jPoutanen fc Svensson 
1996h . 



Numerical calculations of the expected polarization properties of Compton-scattered 
radiation emitted by such hot spots on neutron stars have been done, including special 
relativistic boosting and aberration, and general relativistic light bending in Schwarzschild 
space-time. For slowly-rotating neutron stars, the polarization vector lies in the plane defined 
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by the normal to the hot spot and hne of sight, as expected, since these are the basic 
directions in a non-magnetic system. In fact, the polarization angle can be described by 
a rotati ng- vector type model (see S ecJSl) which includes the special and general relativistic 
effects (IViironen &: Poutanenll2004l ). Model pulse-phase profiles of J, Pl, and x from such 
calculations are shown in Fig{TOl The idea in this case is also to have an atlas of these 
profiles for a grid of values of the essential parameters like spot colatitude (referred to the 
rotation axis of the neutron star) and observer angle, to be compared with GEMS data on 
AMXPs. Such pulse phase polarimetry will be an excellent diagnostic probe of the emission 
geometry in these systems. 



6. Low-Mass X-ray Binaries (LMXBs) 

Bright, accretion-powered galactic-bulge LMXBs are among the brightest X-ray binaries 
known, and are widely believed to be the nurseries producing AMXPs by recycling. Their 
persistent X-ray emission does not show periodicity at any possible rotation period relevant 
for neutron stars, although they are thought to harbor fast-rotating neutron stars with 
millisecond periods. Various explanations have been offered for this, e.g., (a) that the emitted 
pulses are "washed out" by the dense plasma surrounding the neutron star, or (b) that their 
rotation and magnetic axes have become aligned by accretion torques, so that they actually 
do not pulse. Thus, the kind of pulse phase polarimetry we have discussed so far in this 
paper is not possible for them. 



6.1. Burst Oscillations 



However, many LMXBs undergo X-ray bursts, and coherent pulsations at millisecond 
periods have been detected during these bursts in about twenty of them so far. Identifying 
this periodicity with that of the neutron star's rotation, we can again use the above ideas 
of pulse phase polarimetry as a diagnostic probe of the emission geometry in these systems. 
These bursts are thought to be due to the ignition of thermonuclear reactions in the accreted 
material upon reaching a critical temperature and density, deep in the atmosphere of the 
neutron star where the optical depth is extremely large, so that the semi-infinite slab limit 
discussed in Secj5] applies. Furthermore, at sub-Eddington luminosities, the burst process is 
not expected to "lift" the atmosphere to any significant extent, so that the static slab limit 
still applies, and straightforward calculations can and have been done with the aid of the 
Chandrasekhar-Sobolev formalism mentioned above. 
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Results of such calculations (jViironen fc Poutaneru |2004J ) are shown in FigiTTj again 
detaihng the model pulse-phase profiles of /, Pl, and x- The procedure here is again to have 
an atlas of these profiles for a grid of values of the essential parameters like the colatitude of 
the of the spot where the burst occurs, and the observer angle, to be compared with GEMS 
data for diagnostics of the burst emission region. 



6.2. Quasiperiodic Oscillations (QPOs) 

A major property of bright LMXBs in the time or frequency domain is the quasiperiodic 
oscillation (QPO) exhibited by them. QPOs appear as relatively wide, but clearly discernible, 
peaks in the power spectra of LMXBs, as opposed to the extremely sharp peaks corresponding 
to periodic pulses that appear in the power spectra of APPs. QPOs occur at both low 
frequencies {fqpo ~ 6 — 50 Hz) and high "kilohertz" frequencies iyqpo ~ 400 — 1 000 Hz), and 



they have also been detected recently at intermediate "hectohertz" frequencies (Ivan der Klis 



20061 ) . The low-frequency (LF) QP Os are particularly strong, and so were the first to be 



discovered (Ivan der Klis et a/.lll985l ). The study of QPOs and their relations with LMXB 
spectral states has become a most valuable diagnostic tool for probing the dynamics of inner 
accretion disks and the interactions between neutron stars and accretion disks. Studies of the 
polarization properties of QPO sources have a great potential for adding a new dimension 
to this probe. 

Instead of the pulse phase polarimetry methods available for periodic sources, the appro- 
priate diagnostic approach for QPO sources would be through studies of the cross-correlations 
between intensity and polarization parameters, and the associated cross-spectral analysis, 
which, in a sense, are the natural extensions of pulse-p hase polarimetry to qu asiperiodic 



sources. Schemes for such analysis are being constructed (iGhosh fc Swankll2013l ). The ana- 
log of the atlases for periodic sources referred to in the earlier sections of this paper would, 
in this context, be the atlases of cross-spectra for various QPO models, for a grid of values 
of the essential parameters of each model. Comparison of such atlases with the observed 
cross-spectra obtained from GEMS data would discriminate between models and constrain 
model parameters. As LF QPOs are particularly strong, these would be the ideal first testing 
grounds for this scheme. 
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7. Discussion 

We have summarized here the principal polarization properties of the known classes of 
neutron-star X-ray sources, as expected from model calculational schemes summarized in 
the earlier sections. Comparison of the atlases of these expected polarization properties with 
GEMS data will address some key scientific questions on X-ray emission from neutron stars, 
which we now summarize. 

For APPs, the most immediate returns will be on the pencil/fan beam issue. Whereas 
beam-shape diagnostics from pulse-profile studies have been indirect and sometimes am- 
biguous, pulse phase polarimetry provides a diagnostic, viz., the phase-relation between the 
intensity and the polarization degree or angle, which is simple, qualitatively different for pen- 
cil and fan beams, and requires only a coarse binning in pulse phase to achieve clear results. 
However, the most fundamental discovery for APPs would be the discovery of the signature 
of vacuum resonance. As explained in Sec|21 the above phase- relation for a given beam-shape 
changes sign at E^. When E^j lies in the GEMS band, even a very coarse energy-binning will 
be able to find this signature. 

For RPPs, the major scientific issue to be addressed is the nature and site of the X- 
ray emission, as described by the various outer-magnetosphere models summarized in SecEl 
Together with the discrimination between these models through their polarization signature, 
it may also be possible to obtain valuable diagnostics of the magnetic field geometry, e.g., 
whether the field is more similar to a retarded vacuum dipole (Deutsch field) configuration, 
or to a force-free configuration. 

For magnetars, the key issue is again the vacuum resonance point, if it is indeed within 
the GEMS band for the observed magnetars. In contrast to the phase-relation which changes 
sign at Ey for APPs, it does no t do so for the p rimary radiation emitted from the neutron- 



star atmosphere in magnetars (ILai fc Hoi 120031 ). While the impact of propagation effects 



on this diagnostic feature need to be understood fully, all indications at the present time 
are that they do not change it. This is thus a valuable signature of the super strong field 
regime of magnetars, whose demonstration would impact the fundamental physics (QED) of 
vacuum resonance. 

For AMXPs, the key subject amenable to polarizaton studies is the character of non- 
magnetic Compton scattering in neutron-star systems, for comparison with that of similar 
scattering from accretion disks is black-hole systems, particularly at those radii in the latter 
systems where GR effects are mild. This bears on the important issue of the possible impact 
of the remnant magnetic field in AMXPs on the scattering and polarization signature. This 
magnetic field is still invoked for envisaging magnetic chanelling and polar hot-spots of X- 



- 16 - 



ray emission, and indeed one of the model parameters obtained by comparing the atlases 
with GEMS data would be the angle between the magnetic and rotation axes. Any residual 
signature of this magnetic field on X-ray polarization is therefore an intriguing possibilty. 

For LMXB burst oscillations, a key question is the site of thermonuclear ignition and 
X-ray bursts on the neutron-star surface. Recent observations suggest various possible sites 
between the rotational equator and pole for the onset of ignition. Polarization diagnostics 
would give us the angular position of this site in relation to the rotation axis, thus clarifying 
how an X-ray burst develops as the burning front spreads on the rotating star. Once again, 
the role of the residual magnetic field in LMXBs remains an interesting issue here. 

Finally, a key target of polarization studies in LMXBs with QPOs woud be to discrimi- 
nate between various QPO models, particularly those aspects of the models which deal with 
the interaction of the accretion disks in these systems with the remnant magnetic field in 
neutron stars in LMXBs. The role of rotating clumps or blobs of plasma in quasiperiodic 
X-ray emission has been studied intensively in this context, and polarization diagnostics will 
add a new dimension to these studies. 



8. Epilogue 

GEMS was proposed in response to a NASA SMEX announcement of opportunity in 
December 2008. It was selected for phase A development in 2009 and down-selected for phase 
B in 2010. A technically successful Preliminary Design Review was held in Feb 2012. NASA 
Science Mission Direectorate (SMD) indicated their intention to non-confirm (or cancel) in 
May 2012; the SMD decision was based on concerns that the eventual cost would be too 
high. 
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Fig. 1. — Period - Magnetic Field (P - B) diagram of neutron stars, showing regions con- 
taining accretion powered pulsars (APPs, examples denoted by asterisks), rotation powered 
pulsars (RPPs, examples denoted by dots, those in binaries encircled), and magnetars. Re- 
cycled RPPs include both millisecond radio pulsars (MSPs) and accreting millisecond X-ray 
pulsars (AMXPs). Several well-known binary systems and single recycled pulsars are labeled. 
After Ghosh 2007. 
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Fig. 2. — Allowed region, between the two solid lines, in the B — L plane, where the vacuum 
resonance energy falls in the GEMS X-ray band (see text). Only the schematic form is 
shown, with a rough scaling given by standard accretion theory. The upper line corresponds 
to the lower end of the GEMS band, and the lower line to the upper end. 
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Fig. 3. — Atlas of pulse profiles (top panel), polarization amounts (middle panel), and 
position angles (bottom panel) for APPs with pencil beams. In each panel, the five columns 
correspond to different geometries, labeled by the value of ii/i2, where ii — observer angle 
referred to the rotation axis, i2 = inclination angle between magnetic and rotation axes, 
both in degrees. The eight curves in each panel and column correspond to different photon 
energies, labels 1 through 8 referring to energies of 1.6, 3.8, 9.0, 18.4, 29.1, 38.4, 51.7 and 
84.7 keV, respectively. From Meszaros et al. 1988. 
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Fig. 4.— 



Same as FigJSl but for fan beams. From Meszaros et al. 1988. 




Fig. 5. — Pulse profiles, polarization amounts, and position angles for the Crab pulsar, 
showing optical observations from Slowikowska et al. (2009) compared with polar cap, slot 
gap/two-pole caustic, and outer gap models (see text) from Dyks et al. (2004) and the 
striped wind model from Petri & Kirk (2005). The latter model is calculated for two Lorentz 
factors, 20 and 50, shown as solid and dashed lines, respectively, a — inclination angle 
between rotation and magnetic axes, C, — observer angle with respect to the rotation axis. 
Prom Slowikowska et al. (2009). 
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Fig. 6. — Atlas of pulse profiles (black), polarization amounts (blue), and position angles 
(red) for RPPs with a vacuum (Deutscli) magnetic field configuration in the two-pole caustic 
(TPC) model (see text). Angles a and ( as defined in caption of FigJSl From Harding 2012. 
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Fig. 7. — Same as Fig|6l but for the outer gap (OG) model (see text). From Harding 2012. 
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Fig. 8. — Pulse profiles and polarization amounts for magnetar emission from neutron-star 
atmospheres. Line style coded by photon energy: IkeV (dashed), 2 keV (long-dashed), 5 
keV (solid). 7 = observer angle referred to the rotation axis, /3 = inclination angle between 
rotation and magnetic axes. From Lai and Ho 2003. 
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Fig. 9. — Pulse profiles, polarization amounts and position angles for the final magnetar emis- 
sion after including propagation effects in the magnetosphere. Results shown for two energy 
bands: 2-4 keV and 4-12 keV. Line style/color coded by geometry, i.e., the combination (in- 
clination angle, observer angle), as follows. (45,70): solid/black, (70, 45): dotted/red), (60, 
70) : dashed/blue, and (90, 90) : dot-dashed/green. All angles in degrees. From Fernandez 
and Davis 2011. 
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Fig. 10. — Pulse profiles, polarization amounts and position angles for AMXPs. 9 = inclina- 
tion angle between rotation and magnetic axes, i = observer angle referred to the rotation 
axis. From Poutanen 2010. 
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Fig. 11. — Pulse profiles, polarization amounts and position angles for burst oscillations 
from LMXBs. 9 — colatitude of the ignition site and i — observer angle, both referred to 
the rotation axis, i is taken as 60 degrees in all cases shown. Prom Poutanen 2010. 



